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1. Introduction 

One of the main purposes of Statistical Physics is to understand the macroscopic 
behavior of microscopic systems. For regular matter, composed of negative (elec- 
trons) and positive (nuclei) charges, this question is highly non trivial because of 
the long range of the Coulomb potential. 

In 1966, Fisher and Ruelle have in [14] raised the important question of the 
stability of many-particle systems at the macroscopic scale. This may be formulated 
by requiring that the energy per particle (or the energy per unit volume) stays 
bounded from below when the number of particles (or the volume \D\ of the sample) 
is increased, 

\D\ > °- 

For many-body systems interacting through Coulomb forces like in ordinary matter, 
the first proof of stability is due to Dyson and Lenard [TTJIH]. The fermionic nature 
of the electrons is then important [TU] ■ A different proof was later found by Lieb 
and Thirring in [27], based on a celebrated inequality which now carries their name. 
We refer the reader to [22, 23 , 26 for a review of results concerning the stability of 
matter. 

The stability of matter as defined by Fisher and Ruelle only shows that the 
system does not collapse when the number of particles grows. A more precise 
requirement is that the energy per particle (or the energy per unit volume) actually 
has a limit when the number of particle (or the volume |D|) goes to infinity 

lim m = f . 

\D\^oo \D\ 

For short-range interactions, this was already done by Ruelle [5Q] [S3 122] and 
Fisher [15] . The first proof of a theorem of this form for Coulomb systems is due to 
Lieb and Lebowitz in [24] . In this latter work rotational invariance plays a crucial 
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role. For quantum crystals, in which the nuclei are classical particles clamped on a 
lattice (a system which is obviously not rotationally invariant), the first proof goes 
back to Fefferman [TJ]. The main challenge of all these works was to find an ade- 
quate way to prove the existence of screening, the fact that matter spontaneously 
organizes in a locally neutral way. Screening is at the origin of a faster decay of the 
interactions between the particles and it is the main explanation for the existence 
of such systems at the macroscopic scale. The importance of screening was already 
stressed in a fundamental paper of Onsager . 

In two recent papers [TBI HZ], Hainzl, Solovej and the second author of this article 
have proposed a new method for proving the existence of the thermodynamic limit 
for quantum Coulomb system. This method is based on an inequality quantifying 
screening due to Graf and Schenker [15] . and which was itself inspired of earlier 
works by Conlon, Lieb and Yau [8] . The purpose of the present work is to extend 
the results of |16[ 117) to stochastic systems in which the electrons are quantum and 
the nuclei are random classical particles. 

It has been known for a long time that the presence of disorder can strongly 
influence the behavior of a quantum system. The most famous example is of course 
the so-called Anderson localization [2 of particles under weak disorder. On the 
mathematical side, lots of works have been devoted to the study of noninteracting 
disordered quantum systems, for instance described by random Schrodinger opera- 
tors (for an introduction to these results, see, e.g., [Ill [19]). To our knowledge, the 
mathematical literature on interacting disordered many-body systems is quite lim- 
ited, in spite of the increasing physical interest devoted to such systems (3] [9] [33] . 
In recent works (6[ [1] localization bounds were derived for systems of a finite num- 
ber of particles with short range interactions. Some authors considered nonlinear 
random models describing condensed bosonic systems, mainly in Gross-Pitaevskii 
theory (see, e.g., [20] and the references therein). 

In a recent paper [35] . Veniaminov has initiated the mathematical study of the 
thermodynamic limit of random many-body quantum systems with short range 
interactions, following the approach of Ruelle and Fisher. His work does not cover 
ordinary matter made of Coulomb charges, however. Large stochastic Coulomb 
systems were considered before by Le Bris, Lions and the first author of this paper 
in [5]. There the electrons are only described by Thomas-Fermi- type theories and, 
in this case, it is possible to identify the thermodynamic limit / exactly. 

Following the method of [JJ)J [T7] , we are able to deal with quantum electrons 
satisfying the full many-body Schrodinger equation, in the Coulomb field of a ran- 
dom distribution of pointwise classical nuclei. Similarly as in [5], we typically think 
of a perfect (periodic) lattice of nuclei whose location and charges are perturbed 
randomly. For the existence of the limit, we have to assume that this random- 
ness has some translation invariance. This is reflected in the assumption that the 
distribution of nuclei is stationary and ergodic, as we explain below. 

In the next section we properly define our model and we state our main theorem. 
In short, it says that the thermodynamic limit exists and is deterministic, that is, 
independent of the randomness uj: 



(1) 
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More precisely, the limit ([TJ holds in some LP space with respect to the randomness 

(2) lim E 

Almost-sure convergence is expected as well, but not proved in this paper. In 
general, the convergence cannot be uniform with respect to w. In Section [3] we 
will show that if we attach independent harmonic oscillators to the nuclei of a 
cubic lattice, and make them vibrate randomly according to the associated Gibbs 
measure, then 

E(J"(-,L>) 3 ) = +00 

for any D large enough. Therefore, @ cannot hold for p 3 in general. 

For the proof of ([2]) , we will rely heavily on the machinery introduced in [16j EI] 
for deterministic systems. Some parts of the proof which are similar to those of 
PH Q7] will only be sketched. 

Acknowledgement. The authors acknowledge support from the French Ministry 
of Research (ANR-10-BLAN-0101). M.L. acknowledges support from the Euro- 
pean Research Council under the European Community's Seventh Framework Pro- 
gramme (FP7/2007-2013 Grant Agreement MNIQS 258023). 

2. Main result 

2.1. Random distribution of the nuclei. We consider a fixed discrete subgroup 
££ of R 3 , with bounded fundamental domain W . The whole space R 3 is the disjoint 
union of the sets W + j for j e Jzf. We typically think of Jz? = Z 3 with W = 
[—1/2, 1/2) 3 the semi-open unit cube. 

Our nuclei are placed in random locations in the whole space R 3 . We assume 
that the probability distribution of their positions and charges in a given domain D 
is the same when D is translated by a vector k £ . The appropriate mathematical 
notion is that of stationarity which we now recall. 

Let (fi, P) be a probability space. We assume that the discrete group Jzf acts 
on Q and we denote this action by for k G jSf. In the whole paper, the group 
action is supposed to be measure preserving, 

(3) VfcGJz*, P(r fe T)=P(T), 
and ergodic 

(4) (r k T = T,VkeJ?) => P(T) e {0, 1}. 

We now follow the notation of [16] and describe our nuclei by a countable set 
K = {(R, z)} C R 3 x [Z_, Z] with < Z_ < Z. We always make the assumption that 
the nuclei have a highest possible charge Z. Also, their charge cannot be smaller 
than Z_. In reality the charges of the nuclei are integers and they are smaller than 
118. In our random setting, the set K, is random, that is, it depends onwell: 

K(u) = {(i^(w),z,-(w)), R^u) € R 3 , Zj {u) € [Z,Z], j e N}. 

The specific choice of a numbering of these nuclei by the index j has no real im- 
portance. Our main assumption is that the sets K(uS) are stationary with respect 
to the action of on ft, in the sense that 

(5) VfcGJ^, £(t*w) =K(u) - k := {{R- k,z) : (R,z) G JC(u)}. 



\D\ 



-f 



= 0, p>l. 
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A similar setting is used in [5]. Throughout the paper we will always make the 
hypothesis that the number of nuclei in any given set D is bounded almost surely, 
and that two nuclei can never be on top of each other. By stationarity, this means 
that the random variable 



(6) 



X (u}) ;=#{(R,z) e/C(w) : ReW} 



is almost surely finite. It will be convenient to introduce the distance to the nearest 
neighbor of each nuclei in /C(cj), which is the random variable 

R'\. 



(7) 



V(R,z) e K(u), 



inf \R- 

(R' ,z')ElC(u) 
R'^R 




Our assumption that all the nuclei are different means that inf{5R j2 (cL>) : (i?, z) € 
JC(lo) n W} > almost surely. Since the number of nuclei is locally finite almost 
surely, this is the same as asking that the random variable 



(8) 



is almost surely finite. In the following we will write for simplicity (R, z) £ K,(u>)P\D 
to say that (R, z) G K,(ui) and R € D. 

Example 1 (The i.i.d. case). The simplest example to keep in mind is that of 
nuclei on a lattice which are perturbed by independent and identically distributed 
(i.i.d.) random variables (see Figure^). We explain this for a cubic cristal Jz? = 
T? with exactly one nucleus per unit cell. The extension to the general setting is 
straightforward. 

Let us fix a probability space (fio, ^,Po) and consider the product space 



n = (o ) 



ST = a{%) 



= (Po) § 



We choose for the action of Jzf on ft the shift tj, [(w^-gza] = (^j+k)jez 3 - It * s 
known to be ergodic. Consider then r : Qq — > R 3 and z : flo — > [Z_, Z] two fixed 
random variables. The families of i.i.d. random variables 

r j{u) := r(ojj), Zj(oj) := z(u>j) 

are stationary in the sense that rjirki^j) = r{uij + k) = rj + k{u) (and a similar prop- 
erty for Zj). Finally, we let 



(9) 



JC(w) = {(j + r j (u),z j (Lj)) : jeZ 3 }. 



It is obvious that K, is stationary in the sense of ([5]) . We typically think of a 
Gaussian random variable r whose law is given by 

1 „M 2 



(10) 



v{x) = 



(2wa) 3 / 2 

and which corresponds to independent harmonic vibrations of the nuclei. 
The number of nuclei in W is given by 
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We have by stationarity 

P(j + rj eW) = P(r eW-j)= P (r € W - j) 

and therefore 

E(x ) = p 0' + oe^) = E p °( r e w - i) = 1 

Hence Xq G L 1 (f2) is finite almost- surely. The average number of nuclei per unit 
cell is the same as in the deterministic case, it is independent of v. It is possible 
to give conditions on the random displacement r which ensure that Xi is finite as 
well (see Section\3\). 

In this example we have assumed that no nucleus is ever removed from the system. 
The opposite case can be handled by allowing z(uj) G {0} U [Z_, Z] and adding the 
assumption that zj(uj) ^ in the definition of K{w). Then E(Xq) = Fo(z 7^ 0). 




Figure 1 . The case of nuclei on a lattice whose charges and po- 
sitions are perturbed by i.i.d. random variables. A Gaussian dis- 
placement ([TO]) corresponds to having the nuclei attached to har- 
monic oscillators vibrating randomly and independently, according 
to their Gibbs measure. 



As we will see later, in the general case the conditions that Xq and X\ are finite 
almost-surely are not at all enough to prove the existence of the thermodynamic 
limit. But, before writing a more precise condition, we turn to the description of 
the quantum electrons. 

2.2. The grand canonical free energy of the electrons. Here we mainly fol- 
low [T7] . Let D be a bounded open subset of R 3 . The Hamiltonian for N electrons 
in D is the random self-adjoint operator 

N N 

(11) H(u,D,N) = Y(-A) x , v " 



' ** ' ^—^ \x„ — R\ 

n=l n=l (R,z)eK(u)nD 1 1 



\x n -x m \ 2 ^ \R-R' 
{R' ,z')eK{uj)r\D 

R^R' 
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Here —A is the Laplacian with Dirichlet boundary conditions in D. The operator 
H{uj,D,N) acts on the fermionic space L 2 {D) consisting of square-integrable 
functions ^(x\, xn) which are antisymmetric with respect to exchanges of the 
Xj. The form domain of H(ui,D,N) is the Sobolev space A, H^j(D). Since there 
is a finite number of nuclei in D, which are all distinct, for almost every to G f2, 
the operator H(uj, D, N) is self-adjoint on H 2 (D) n Hg(D) and bounded from 
below. For simplicity we have chosen units in which h ~ 1 and the mass m and 
charge of the electrons are m = 1/2, e = 1. We have also neglected the spin for 
convenience. The terms in H(uj, D, N) respectively account for the kinetic energy 
of the electrons, the nuclei/electrons attraction, the electrons/electrons repulsion, 
and the nuclei/nuclei repulsion. 

At zero temperature, the ground state energy for N electrons in D is the random 
variable 



(12) F q (uj,D,N) = mia(H(uj,D,N)) 



inf 

*eAf Hl{D) 
f D N l*l 2 =i 



(<£, H{co,D,N)^). 



Since H(uj,D,N) is bounded from below, JFo(ui, D, N) is a well-defined random 
variable. As we will explain, the stability of matter tells us that there is a lower 
bound on J-o(uj, D, N) which is independent of u> and N (see Theorem[2]below). The 
stationarity of the distribution of the nuclei implies a certain stationarity property 
of Jo which reads ^(t/cW, D, N) = F^lo^D — k,N). Indeed, the Hamiltonians 
H(rkUj, D, N) and H(u, D — fc, N) are isometric and therefore they have the same 
spectrum. 

Like in [121 117] . we work in the grand canonical ensemble, that is, we optimize 
over the number of electrons, instead of imposing the neutrality of the system. The 
ground state energy in D for a realization of the distribution of the nuclei is then 
defined as 



(13) 



To(u,D) := w£ o ? (u,D,N). 



This random variable is stationary in the sense that 

(14) F (TkU,D)=F (u,D-k) 

for any chosen domain D, all k € JSf and almost all uj G ft. 

At positive temperature T > with chemical potential (i6l, the grand canon- 
ical free energy is defined by the formula 



(15) 



7t, p (w,D) = -T log ^Tr A? 



N r,2fB3-| e 



-(h(u,D,N)-i*n) /T 



,N>t) 



and it satisfies a similar stationarity property as J-q- 

There is a useful variational formula for Tt,^ in the fermionic Fock space & 
C © 0jv^i Ai^ L 2 (D). Introducing the operators 

H(w,D) :=0© H(u,D,N) and A^:=Oe0JV, 



N>1 



N>1 
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we then have 

(16) FtAu> d )= inf Tr&((W(u},D)-ijM)T + TTlogT 

Tr^r=i 

Here T is the density matrix of a mixed quantum state in Fock space. See |17j for 
more details. 

2.3. Existence of the thermodynamic limit. Before stating our main result, 
we quote the following important lower bound on Tt^i which is nothing else but 
the stability of matter. 

Theorem 2 (Stability of matter). There exists a constant C depending only on 
the highest nuclear charge Z such that 



(17) 



F Ttli (u>, D) > -C(l + T 5/2 + a4 /2 ) \D\ 



for all T 0, fi £ R and almost all 

Here fi + — max(0, /i) denotes the positive part of fj,. We see that J-t, ai (w, D) can 
be very positive (depending on the number and on the positions of the nuclei), but 
it can never be too negative. Theorem [2] was proved for the first time by Dyson 
and Lenard in [11[ 121] and it was later revisited by Lieb and Thirring in |27j . For 
a recent proof, see [26] and [17j Theorem 3]. The only important property of the 
nuclei used to get this bound is the fact that their charge is uniformly bounded by 
Z. The stationarity and the fact that the number of nuclei is locally bounded are 
not used to get the lower bound (fT7| . 

If the charges of the nuclei are not random but they are all equal to the same 
charge Z, there is a better lower bound on Jr,^(w, D): 



^ 5r, z 



(18) Tt^uj, D) ^ -C(l + T 5 / 2 + ft 2 ) \D\ + — 

(fi,2)eK(u)nD 

where we recall that <5r jZ is the distance of (R,z) to the nearest nucleus in )C(lo). 
This bound can be obtained by using an inequality due to Lieb and Yau [25] (gen- 
eralizing another one of Baxter [4]), see [26] and [TT] Theorem. 3]. In the thermo- 
dynamic limit, we hope to prove that J-t,/h(w, D) ~ C\D\. A necessary condition is 
at least that 

e( 

for D smooth and large enough. By stationarity, this is equivalent to 

(19) E{X 1 )=eI J2 J- 

\(R,z)eK(u>)nw R,r 

We recall that Xi is defined before in ([5J. This inequality says that inverse of 
the smallest distance between nuclei in W has to be summable in average and this 
implies that the number Xq(uj) of nuclei in W (defined in ([6])) has to be in L 1 (fi) 
as well (see Lemma [TH below 1 ) . All this is already much stronger than saying that 
the number of nuclei is finite and that the nuclei are distinct almost surely, as we 
have done before. 



< oo. 
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We are able to prove the existence of the thermodynamic limit under a similar 
but stronger assumption than (p~9|) . see ([20)) . Before stating our main result, we 
however need to introduce a notion of regular domains, following [161 117) . 

Definition 3 (Regular domains [El [17]). 

• (Sets with regular boundary) Let a > be a real number. We say that a domain 
D C K 3 has an a-regular boundary in the sense of Fisher if 

Vie[0,l/a), | a; £ M 3 | d(x,dD) |L>| 1/3 tj ^\D\at, 
where 3D = D \ D is the boundary of D. 

• (Cone property) Let e > be a real number. We say that a set Ac I 3 has the 



e-cone property if for any x £ A there is a unit vector v x £ 



such that 



{y £ K 3 | (x - y) ■ v x > (1 - e 2 )\x - y\, \x - y\ < e} C A. 

• We introduce the set 1Z. a ,e of all bounded open subsets D C R 3 which have an 
a-regular boundary and such that both D and R 3 \ D have the e-cone property. 

Our main result is the following. 

Theorem 4 (Existence of the thermodynamic limit for random quantum Coulomb 
systems). We assume that the distribution of nuclei K,(lj) is stationary in the sense 
of (|5|) and that it satisfies 



(20) 



ll*i I 



LP (SI) 



E 

(R,z)£K{ui)nW 



1 



< oo 



LP (SI) 



for some 2 ^ p ^ oo. Then there exists a ( deterministic) function f(T, /i) such that, 
for any sequence (D n ) C TZ a ,e with a, e > 0, \D n \ — > oo and diam(_D„)|Z? n | -1 / 3 ^ c, 
we have 



(21) 



lim E 

n— >oo 



\D„ 







for q — 1 if p — 2, and for all 1 q < p/2 if p > 2. In particular, we have for a 
subsequence 



(22) 
almost surely. 



lim 



\D T . 



We believe that almost-sure convergence as in (|22[) holds for the whole sequence 
D n , provided that it does not escape too fast to infinity. A simple condition is that 
D n C -B c i£) 11/3 where B r denotes the ball of radius r, centered at £ M 3 . A condi- 
tion of this type is needed for Birkhoff's almost-sure ergodic theorem (Theorem [22] 
below). Actually, under this additional assumption, we are able to prove that 

J"t, m (w,Ai.) 



(23) 



lim inf ■ 



\D„ 



see Remarks 1241 and 1261 below. However, proving that (|23[) holds with a limsup may 
require more involved tools from the theory of probability. 

In the next section we will give simple examples of distributions of nuclei satisfy- 
ing the condition ([20]) . Recall that when all the charges are equal to Z, a condition 
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of the same form as (|20|) is necessary, with the power p = 1. Our assumption (|20|) 
with p ^ 2 on the distribution of nuclei is far from optimal. The power 2 (which is 
due to several rough kinetic energy estimates) is probably an artefact of our proof. 
We have not tried to improve the condition ([20]) too much. 

It is possible to prove that, in the thermodynamic limit, the system wants to be 
neutral in average. This can then be used to show that the chemical potential (i 
does not play any special role here: f(T,fi) is just linear with respect to /z. 

Corollary 5 (Asymptotic neutrality and form of f(T, fi)). Under the same as- 
sumptions as in Theorem^ let 



Eat -(H(u,D n ,N)-uN)/T 



N Tilt (u,D n ) 



V" Tr „ -(H(u,D n ,N)-nN)/T 
AT>0 

be the total average number of electrons in D n . Then we have 



(24) lim 

n— )-oo 

where 







(25) Z av := — E\ £ : 

1 1 \(fi,z)effnt(i)) 

is the average nuclear charge per unit cell. Moreover, we have 

(26) f(T,n) = f(T,Q)-nZ m 
and T /(T, 0) is concave. 

It is easy to see that Z av is finite under our assumption (|2"0)l on X\ , see Lemma[T4l 
below. Corollary [5] is proved later in Section[5] It easily follows from the stability of 
matter (Theorem^ and the upper bound E(.Fr,/j(', ^ C| Z? T1 1 which is proved 
in Lemma [T2] below. 

Remark 6. In this paper we have considered a discrete group acting on R 3 and on 
the probability space Q, in an ergodic fashion. Our method of proof can be applied 
to deal with a continuous group like K 3 . For instance, the thermodynamic limit 
exists if the nuclei are distributed using a Poisson process, since the corresponding 
random variable X\ satisfies (I20p in this case. 



3. The case of i.i.d. perturbations of a nuclear lattice 

In this section we come back to the special (but instructive) case of i.i.d per- 
turbations of nuclei on a lattice, as introduced in Example [1] In this setting we 
derive conditions under which the assumption (|20[) on X\ in Theorem U is satisfied. 
For simplicity we consider a cubic lattice Z 3 with only one nucleus per unit cell 
W — [—1/2, 1/2) 3 . The extension to a more general situation is straightforward. 

As in Example [U we assume that Q = (fl ) z3 , P = (Po)® 2 ' and that 

(27) /CM = {{j+ r» , Zj (cj)) : j el?). 
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with Vj(uj) = r(u)j) and Zj(u>) = z(u>j). Here r : £Iq — > M 3 and z : fto — > [Z_, Z] are 
two fixed random variables. Let us denote by v the law of r. We will give conditions 
on v which ensure that X\ G L P (VL) for some p ^ 2, as required in Theorem^ The 
simple Gaussian case (|10[) is covered by the next two results. 

We recall that the random variables Xq and Xi are defined by ([6]) and ([8]), 
respectively, that is, 

X (u) = # (/C(w) n W) = # {j e Z 3 : j + r^w) G M^} = E 1(7 + 'jM e 

jez 3 

and 

= ^ iCj + ^He^) ^ w . th = b . + _ k _ _ 

Recall that P(j + r } G W) = P (r G W - j) = u(W - j), hence 

E(X ) = ]T P(j + r,- G W) = "( W - 3) = v(^ 3 ) = 1 

j£Z 3 jSZ 3 

and in particular Xq G L 1 (17). The following elementary proposition deals with the 
intcgrability of higher powers of Xq ■ 

Proposition 7 (Integrability of Xq in the i.i.d. case). We assume that JC is of the 
form (|27p. Then 



(28) ll*o|| iP( o) < $>(W-j) 



i/p 

jez 3 



T/ims Xo G -L p (f2) w/iera i/ie rzg/ii side is finite. If the support of is is not compact, 
then X Q <£ L°°(f2). 



Proof. Using that X = X^- e z 3 1 + r J e ^0 » we obtain 

i^oi LP( o) < E e w )IL,(n) = E p t?+^ e = E <W-j) 1/v . 

jez 3 jez 3 jez 3 

This bound does not really use the independence of the variables . We now prove 
that Xq ^ L°°(0) when the support of 1/ is not compact. This means that there 
exists an infinite sequence (J n ) C 1? such that 

Vn, P(r in £W — j n ) = P (r eW- j n ) = v(W - j n ) > 0. 

Let TV be a positive integer. It is clear that if for all 1 ^ n ^ N we have j n +rj n G W, 
then Xo N. Hence, we have 

(N \ N 

n On + r iB e = n p o» + ^ e 
n = l / n=l 

Since V(j n + r Jn e W) > for all n, we deduce that P(X ^ N) > for any X > 0. 
This implies that X g L°°(0). □ 

The following is now a simple application of the previous proposition. 

Example 8 (Gaussian perturbations I). Assume that v{x) — (27rcr)~ 3 / 2 e~' :r ' 2 ' / ( 2cr ) 
is a Gaussian distribution. Then Xq G L p (Q) for all 1 ^ p < oo 6m< Xq ^ iy°°(r2). 
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We next turn to the study of X\ . The simpler situation is when the nuclei never 
escape from their cell and stay at a finite distance to the other ones. 

Lemma 9 (Small perturbations). If v has its support inside W , then both Xq and 
Xi are m L°°(0). 

Proof. Under the assumption, it is obvious that there is always exactly one nucleus 
per unit cell, and that it is at a finite distance r\ to any other nucleus. Thus Xq = 1 
and Xi^ I/77. □ 

Corollary 10 (Thermodynamic limit for small i.i.d. perturbations of the nuclei). 

We assume that K, is of the form (|27| and that v has its support inside W . Then 
the thermodynamic limit in Theorem^ 

(29) lim E 

n— too 

is valid for any 1 ^ q < 00 . 

A more interesting situation is covered in the following 

Proposition 11 (Integrability of X\ in the i.i.d. case). We assume that JC is of 
the form (|27[) and that v satisfies 

(30) E I 

for some rj > and some 1 ^ p < 3. Then we have X\ £ L P (Q). 

If moreover there exist a ball B K (v) C W with radius n > and center v £ W , 
and i 7^ j £ Z 3 such that 

(31) Vx £ (B K (v) - i) U (B K (v) - j), v{x) ^k>0, 
then Xt i I 3 (Q). 



The condition (1301) implies that the measure v is actually a bounded function 
outside of the unit cell W (indeed, outside of (1 — 7]/y/3)W), and that it decays fast 
enough at infinity, in a similar fashion as in (|28p . On the other hand, v does not 
have to be absolutely continuous with respect to the Lebesgue measure inside W. 
The condition pip means that, with a positive probability, two nuclei coming from 
different sites i and j will both (independently) visit all of the same ball B K (y). 

Proof. Recall that 

Xi = E 1 ^(i + r i)^ 



therefore 



(32) 


I^ll^(n) < E 


1w0' + ^)t- 


= E 


. . 1 
tw-j{r )-r 




jez 3 


°j 







for p ^ 1. In order to show that the series on the right side is convergent, we first 
estimate P(ro £ W — j n So < e). For this purpose, we point out that 5q < e if 
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and only if there exists k ^ such that |ro — k — r k \ < e. Hence, 

P(r E W - j n S < e) ^^2P(r Q EW - j n |r - fc - r k \ < e) 

= E/ W K»)*U.(.-*-,) 

k7i0 Jw-j Jr3 

k ^ Jw Jb e 

Even if ^ is not necessarily absolutely continuous in the interior of W , we have used 
an integral notation for simplicity. When j = 0, we obtain, for e < r\, 



P(r eW n <5 < e) < V" / i/(x) dx v(y + x - k) dy 
kjt0 Jw Jb e 



< E HIl^+b,,-*) (M3)£ 3 = Ce 3 v{W). 

\k^o J 

In the second line we have used that v(y + x — k) ^ IWl^iw+B^-k) f° r £ < V- Note 
that 

E Hi~(W+.7+l?„) < 00 =^ E I I/ lli«'(lV+j+B„) < 00 

since p ^ 1. When j ^ 0, we isolate the term k = —j and obtain 

-W 

fe/0 "' W ' 

We estimate the first term by (recall that e < rj) 

where we have used v{W — j) ^ IHIl<»(vk-.j)- F° r the second term we write 

/ v{x-j)dxi v{y + x)dy^\v\ LX , w _ j) l dx j v(y + x)dy 
Jw Jb e Jw JB e 

< Ml<~{w-3) v(W + B v )(4ir/3)s 3 

= CM LOO[w _ j) e 3 . 

As a conclusion, we have shown that 

'v{W) forj = 0, 



P(r eW-j n 5 < e) < Ce 



3 



lkli,°o(w-j) for.?' 7^0. 
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We then compute the expectation value of l(ro G W — j) 5q 9 : 



E 



Mro eW-j) 
°o 

t I 1U ^2 

<5o 



E E 



l(r G W-j) 



°0 



< + £ 2« fe P (r <EW — j n (5q < 2-( fc -!) 



2« 

- j) 
2i 



c(i(j = omw) + i(j ? o) ) J2 2^ fe 2- 3 ( fc - 1 ). 



fc>0 



The sum is convergent provided that 1 ^ q < 3. Hence we have shown that 



VI < q < 3, 



l(r e W-i) 



So 



c 



L*(Sl) 



\v{W) 1 l q forj = 0, 
(1 - 2«-3)V ? - {l^^ forj^O. 



Inserting this in (|32[) gives the result under our assumption (j30|) . 

In order to conclude the proof of Proposition [Til we show that X\ $ L 3 if ([31 
holds. For this purpose, we write 



*i- E E E w > E 6 f 

iez 3 jez 3 fcez 3 J iez 3 l 

i+ri£W j+ rj eW k+r k £W i+riGW 



Therefore it is sufficient to prove that S~ 3 t(i + r.t G W) L 1 (51) for some i. We 
choose i,j such that (13"TT) is satisfied. We have in particular P(i + r, G W) = 
i/(W — i) > 0. We compute 



F(5i>e n i + r 4 G W) =P ( {i + G VF} n Q {|i + r 4 - fc- r fe | > e} 

/ \ 



lim 

N—nx 



{i + r t eW}n P) {\i + r t - k - r k \ > e} 



V 



\k\^N 
k^i 



1 
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Using 1 — t ^ e * , wc obtain the bound 

/ \ 

{i + neW}n f] {\i + n - k - r k \ > e} 



V 



\k\^N 



J 



v{yi)dyi Yl / v iVk) dy k t(\i + y t - k - y k \ > e) 

\kKN jR3 

/ v{yi) dy l TT ( 1 - / v{y k + k - y t - i) dyj ) 

Jw-i V Jb c J 



|fc|<JV 
k^i 



< / K2/*) ex P 

JW-i 



B e 



\ 



^2 v(y + k-yi-i) 



i |fe|<iV 



dy 



J 



dyi 



Passing to the limit N — > oo, we have shown that 



"(Si > e n i + n G W) < / v(yi) exp 



/ u (v + k - Vi) \ d v 

Jb - \k^0 J 



dyi. 



When s < k/2 we can use our assumption (|3ip and infer that v(y + k — y^) > k for 
yi G B K / 2 (v) — i, y G B £ (0) and k = j — i. In particular, 

p(6i > e n i + neW) 

< / v(yt)dy t + ^(yi)exp [-Cke 3 ] dy, < (1-Ce 3 ) / v{y l )dy i 

JW\B K (v)-i JB K (v)-i JW-i 

where C > depends on i (a fixed index). We deduce that 

P(Si ^ e n i + n G W) ^ Ce 3 P(i + G W) 

and finally obtain 

'l(i + ri G W). 



E 



<5? 



l(&<e) > CP(t + r 4 G W). 



If l(i + Ti G W)5^~ were in L 1 (J7), the left side would converge to when e — > 0, 
by the dominated convergence theorem. Since the right side is > by our choice of 
i, and independent of e, we deduce that t(i + n G W)5^ 3 ^ This concludes 

the proof of Proposition QTJ □ 

Let us recall the inequality (|18|) which implies that, when z(u) — Z a.s., 

^ ^ liMa^fw.^ + cfl + T^+pf)^!. 
j ■ 

W-jCD 

This clearly shows that Tx^ i^, D) is not in L 3 (f2) when (f5T|) is satisfied. 
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Corollary 12 (Thermodynamic limit for large i.i.d. perturbations of the nuclei). 

We assume that JC is of the form (|27j) and that v satisfies 

( 33 ) EHi»m-j) <0 ° 

for some rj > 0. Then the thermodynamic limit in Theorem^ 

(34) lim E 

n— >oo 

is valid for all 1 ^ q < 3/2. If moreover z(ui) — Z a.s. and (|3ip is satisfied, then 
Jtt,h('i D n ) £ L (fi) /or any n. Ttes (|34[) cannot hold with q — 3. 

Example 13 (Gaussian perturbations II). Assume that v{x) = (2iro-)~ 3 / 2 e~\ x \ 2 / < - 2cr ) 
is a Gaussian distribution. Then we have for all j € Z 3 

^ C e -lilV(^). 

i«(W-i) 

Since (|30l) *s satisfied for all 1 ^ p < 3, we /icwe € L p (f2) /or aZZ 1 ^ p < 3. On 
</ie other hand, the support of v is the whole space, hence (|3ip is obviously verified, 
thus Xi £ L 3 (fl). The thermodynamic limit exists in L q (fl) for all 1 q < 3/2. 
But Jt, m (', D) L 3 (fl) when z(lj) = Z a.s.. 

In the next section we provide the detailed proof of Theorem 0J 



\D„ 



x\ 2 /(2a) 



4. Proof of Theorem [4] 

Our proof follows the technique introduced in [16\ 117) . In |16) , abstract con- 
ditions called (A1)-(A6) ensuring the existence of the thermodynamic limit for a 
functional D H> F(D) were provided. These conditions were verified in [T7] for 
the deterministic crystal as well as some other quantum systems. Our technique of 
proof for the stochastic case can be sketched as follows: 

(i) We start by proving in Lemma 1 151 below that \D n \~ 1 J r T,)j,(<jJ, D n ) is uniformly 
bounded in LP/ 2 (fl). This uses an adequate trial state together with the assump- 
tion (l20l) on the distribution of nuclei. 



(ii) Then, we show that the averaged free energy E (j r T,^(-, D)) satisfies all the 
abstract properties (A1)-(A6) of [TBI [P7]. Hence its thermodynamic limit exists 
and we call f(T,/i) the corresponding limit. Note that E(j r (-,Z3)) is periodic by 
stationarity of J~t,h, hence the formalism of [THJ [T7] is appropriate. This step 
requires some upper bounds in average (in particular the estimate (A4) proved in 
Lemma [TBJ below) in which the assumption (|20p on the distribution of nuclei is 
again used. 

(Hi) We show that 



(35) lim E 



n— >oo 



\D ri 



where [x]— — max(0, — x). This step uses the Graf-Schenker-type inequality (A5) 
(Lemma 1181 below) which is a precise lower bound on Ft,h{w, D) at fixed u>, in 
terms of a tiling of simplices. The limit ()35|) is obtained by a suitable application 
of the ergodic theorem (Theorem [22] below). 
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(iv) A simple argument shows that the convergence of the average E LFt,h(-, D)) 
and (|35p imply the strong convergence in L 1 (J7). By interpolation, the convergence 
in L q (Q.) for 1 ^ q < p/2 follows. 

Actually in our proof we do not consider the original functional Jr^w, D) but, 
like in [17] , we optimize over the charges of the nuclei which are close to the bound- 
ary of D. This provides a modified functional T_ T M (w,_D) to which the previous 
scheme is applied. Only in the end of our proof we come back to the original free 
energy. 

Step 1. Bound on Xq. A preliminary result is the following: 

Lemma 14 (Bounds on X± give bounds on Xq). Assume that the distribution of 
nuclei K{lo) is stationary in the sense of (JSJ) . Then, for any p ^ 1, X\ € L P (Q,) 
implies Xq £ L P (Q). 

Proof of Lemma[T4\ Recall that Xq ^ 0, so we have 

(36) E (\Xq\p) = E (XP) = E (X$l Xo<1 ) + E (X$l Xo>1 ) ■ 

Next, we point out that if Xq > 1, then, for any Re JC H W, we have 8r,z ^ 
diam(H^). Hence, 

Xitxa>i > ^w) lxo>i - 

We insert this estimate into (05]), finding E(|X | P ) ^ 1 + diam(VF) p E(Xf ), which 
concludes the proof. □ 

Step 2. Upper bounds. In this first step we will establish some upper bounds in 
average, that is for |J 7 t, j u( - ; -D)lia(n)- It is f° r these uper bounds that we will need 
the assumption (|20[) which gives estimates (in average) on the number of nuclei per 
unit volume, as well as on the smallest distance between them. The lower bounds 
will on the contrary be almost uniform mu€fi (up to a small error term which is 
easily controlled, see Lemma [T51 and Remark Q1J] below). The proofs of the results 
in this first step are rather technical but the strategy is similar to that used in the 
deterministic case in |16j . 

Our first result will be that, under our assumption (|20[) on the nuclei, the free 
energy is bounded above by a constant times the volume \D\, in average. The 
following is the random equivalent to |16[ Prop. 2] . 

Lemma 15 (Upper bound). Under the hypotheses of Theorem^ we have, for any 
regular domain D € lZ a ,s 

(37) e\t t J,D)\ P/ \c\D\p/ 2 

where we recall that p appears in (|20(1 . and where C depends on a > 0, e > 0, 
T 0, p and /i e R, but not on D € 7Z a ,e- 

In particular, if we denote by T(u>) the electronic density matrix in Fock space & 
of any optimal state for J-T lAt (w, D), we have the bound in average for the electronic 
density pr and the kinetic energy 

C-SN) E j / ! / ,//• '+ Trr( > ! A), r) j <C\D\« 




Tr^(^(-A),r)^ 



for allls^q^ p/2. 
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Proof of Lemma \75l Our proof follows the one of [171 Prop. 2] where some missing 
details can be found. There, the nuclei were all assumed to have a finite distance 
to their nearest neighbor. Our task here is to exhibit the dependence of this upper 
bound in terms of these parameters. 

We only write the proof for /j, = T = 0, the general case being similar [17] . 
and we denote J- = J-q.q. Let D be a regular domain and lj G f2. Recall that we 
have stability of matter which tells us that ^ — C\D\ almost surely. We 

therefore only have to prove an upper bound, which is done by constructing an 
appropriate trial state. For each nucleus (R, z) £ K.(oS) H D, we place in D radial 
electrons of total charge z, in a small ball of radius S' R z /8 where we have defined 
for convenience 

S R,z : = min ( S R,z> £ )- 

Recall that 5r z (cj) is the distance of the nucleus (R, z) to the closest nucleus in 
the system and that e quantifies the cone property of the set D. We want to put 
this ball as close to the nucleus as possible. When the nucleus is at a distance > e 
to the boundary dD, we can simply put the radial electrons on top of the nucleus, 
leading to a vanishing Coulomb potential outside of the support of the electrons, 
by Newton's theorem. When the nucleus is at a distance ^ s to the boundary of 
£>, we use the cone property and place the ball at a distance 5' R z /4 to the nucleus, 
in the small cone which is inside D. This construction is the same as in [17] except 
that our electrons live in small balls depending on 8' R z . In [17 they were all living 
in balls of constant radius e/8. 

To simplify our estimate, we use the notation 

(dD) e := {x e R 3 : d{x, 3D) s? e}. 

Since D has an a-regular boundary, we have |(9-D) £ | ^ ae\D\ 2 / 3 . 

The total energy of our trial state contains several terms. The kinetic energy 
used to squeeze the electrons in their small balls can be estimated by a constant 
times 

V z . 

(Ji,z)EK(u)nD V R - z ' 

Here we get a coefficient z 5 / 3 because of the Pauli principle for the electrons. We 
first have to put z electrons in a ball of radius 1. For this we just fill in the 
first eigenvalues of the Dirichlet Laplacian of the unit ball and then average over 
rotations to make our state radial. Then we scale these electrons to make them fit 
in a ball of radius 5' R z . Here the term z 5 / 3 is not a problem since the charges z 
are uniformly bounded by assumption, z ^ Z. But later this difficulty will pop up 
again. 

The only other term is the interaction between all the charges in (dD) e . The 
interaction between each nucleus and its screening electronic cloud is negative and 
we can discard it for an upper bound. For later purposes, we however note that it 
can be estimated by a constant times 

(39) E 7- 

(R,z)£K(w)n(dD) e R ' z 
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Figure 2. The sets of indices J and dj in the simple example of 
a cubic lattice. 



We are left with the dipole-dipole interactions, which we denote by 

TV ( D tV \ ZZ ' ZZ ' ZZ ' ZZ ' 

P( ' ] = \R-R'\ + \X-X'\ ~ \X-R'\ ~ \X'-R\ 

(even if it also depends on to). Here X and X' are the positions of the electrons 
which are such that \R - X\ < S' R z /A and \R' - X'\ S' R , Z ,/A. When \R - R'\ is 
sufficiently large, this interaction behaves like \R — R'\~ 3 . When \R — R'\ is small, 
we use that the electrons are at a small distance to the nuclei. For instance 

\R X'\ >\R- R'\ \R' X'\ = \R- R'\ - 6' R ,, z ,/4 >\\R-R'\. 

Similarly, \X - X'\ > \R - R'\/2. All in all, we deduce that 

(40) |Dip(i?,i?')| < : rrr. 

y ' \R- R'Kl + lR- R'\ 2 ) 

Our final bound on the energy is therefore of the form 

(41) F{u,D) 

( 



^ C 



5 / 3 ~ 



E mi \2 + E 



(8' V ^ \R- R'\(1 + \R- R'\ 2 ) 

\ RjtR' I 



In order to simplify our reasoning, we now cover D and (dD) e by translations of 
the domain W. This means we write D C Uj^jWj and (dD) £ C Uj^gjWj where 
JCif and dj C if are such that #J ^ C\D\ and #dj s? C|L>| 2 / 3 , by the 
regularity of D (Figure [2]). We also use the notation Wj := W — j. Then, the first 
term of (j4Tj) can be estimated by 

z 5/3 E t/« < ^ £ x '^) < c ipf E 

(R,z)£K(uj)nD K R ' z > JGJ WU JGJ 



by stationarity, and with the definition 
(42) AT» := J2 
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Note that, since (^J -1 < + fez) -1 

1 



E 



(§> )2 

(R,z)GlC(uj)nW y R ' zJ 



E 



E 



K (R,z)eK(uj)nw 
1 + 



J R,z 



(R,z)eK(u)nw (R,z)elC(oj)nW UR < Z 



This shows that X' 2 G L p / 2 (f2), under our assumption ([20|) and using Lemma IT4l 
By the triangular inequality, we get the estimate 



E 



(X' )2 

(R,z)GK(ui)nD y R > z > 



< CPI II^IU* 



Lp/ 2 (n) 



Lp/ 2 (0) 

Note that we implicitly use here that p 2. 

We now claim that the second term of (|4"Tj) can be estimated as follows 



(43) 



E 



(R,z), (fl',y)e/C(a»)n(9£») e 



\R-R'\(1+ \R-R'\ 2 ) 



II LP (O) 



Lp/ 2 (S~2) 

^ c]z»| a / 3 (l-x-ii^^, +iog(i + |z>|) |jr |5 

When R and i?' belong to two adjacent domains Wj and Wk, or to the same domain 
Wj = Wk, we only use that 

1 1/1 1 



\R-R'\ " 2 fo„ 

by definition of 6r iZ - When R and R' belong to two domains Wj and Wk which 
are separated by a finite distance, we use the estimate 



E E 

(R,z)£K.nw j (R',z')eKnw k 



\R-R'\{1 + \R-R'\ 



^ C- 



1 + \j - fc|3 



where we recall that -Xo(w) is the total number of nuclei in the unit cell. The final 
estimate on the second term of (|4T1) is 



(44) 



C Y, X 1 (r J u;)X (TkU J ) + C £ 



X (TjUj) X (T k Uj) 



\j-k\^C 



j,keaj 



+ \j - fc| 3 



We use that 



and obtain 



Xo/i{TjUj) X Q/1 (r k uj) < ^(x 0/1 (TjU}) 2 + X 0/1 (T k u>) 2 ^ 



j.kedj 
|i-fcKC 



Xi(TjU)X (T k U)) 



C\D\ 2 / 3 (\Xi\ 



2 

LP (SI) 



\\Xo\\ 



LP (SI) 



Lp/ 2 (SI) 
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Similarly, 



X (TjUj) X (T k Uj) 



*ZC\D\*/*]og(l + \D\) J2 X (r^) 2 j 



Therefore 



£ 

j,kedj 



X (TjUl) X (T k Ul) 



i + li - fc| 3 



<C|D| 2 / 3 log(l + |D|) IJfolL 



2>(n) 



LJ>/ 2 (fi) 

and (|4"5)) is proved. This concludes the proof of (|57|) . 

The bounds ([55]) follows from the Lieb-Thirring inequality and the stability of 
matter. By (| 1 T[) with a 1/2 in front of the kinetic energy instead of a 1, we see that 
the total energy is bounded from below by 

T T ^(u,D) > I Tr^ ^(-A)ir(w)j - C\D\ 

almost surely. Hence our upper bound on Ft^u, D) yields 



Tr^ (^>A) 4 i»^ 



< C\D\ 



Lp/ 2 (0) 

for a regular domain D. By the Lieb-Thirring inequality we have 



Tr^ ^(-A) t T(o;)^ > C jf > C\ 



£,1-2/3 



PT(u 



5/3 



almost surely, which gives (155)) . This ends the proof of Lemma [15j 



□ 



To simplify some estimates from below that we will derive later, we now introduce 
an auxiliary free energy J T(1 (u),D) obtained by minimizing over the charges z of 
the nuclei (R, z) € K.(ui) D D, which are at a distance ^ 2e from the boundary of D. 
This means we replace the charge z(u>) of each of these nuclei by z 1 and we minimize 
over these z's, under the constraints that ^ z' ^ z(uj). This trick simplifies some 
lower bounds and it was also used in |17) . The idea is to show the existence of the 
thermodynamic limit for T_ T and, only in the end, to prove that this implies the 
result for the original function Ft,h- This is done by using Lemma [TBI below, and 
the fact that 



(45) 



for all D and almost all w. Note that the random variable D) satisfies the 

same stability of matter inequality (|T7|) as Jr ift (w, D), and by (l4"51 it satisfies the 
same upper bound ()37p as Jr i(1 (w, Z?). It is also a stationary function in the sense 
of ([14")) . The following is the equivalent of [TTJ Prop. 4] in the random case. 
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Lemma 16 (Control in average of charge variations at the boundary). Let D' € 
lZ a ' ,e' and D € TZ a ,e be two regular domains such that D' C D and d(dD, dD') C . 
Then we have 

(46) E(T Tjfl (;D)) ^E(Z T J;D'))+C\D\D'\ + C\D\ 13 / 15 

where C > depend on a, a! , e, s' , W , \i and T , but not on D and D' . 

The power 13/15 is not optimal. In the end of the proof we indicate how to 
improve it. 

Proof. Like for the proof of Lemma [T51 above which was based on [T71 Prop. 2], we 
now follow the proof of [TTJ Prop. 4], but we keep track of the smallest distance 
between the nuclei and we use a slightly different argument in the end. To simplify 
our reasoning we assume that e = e' and that T = fj, = 0. We denote T := T . 
The present proof easily carries over to the general case. 

(i) The trial state. For any fixed w G f! we pick the exact trial state T for the vari- 
ational problem >F(-, D'). We then use this trial state to get the upper bound (l4"bT) . 
Recall that for the variational problem >F(-, D') the charges close to the boundary 
of D' arc optimized. In our system the nuclei do not necessarily have these optimal 
charges. In D\ D' we have several nuclei which we have to screen. Like in the proof 
of Lemma [T5l we do this by placing electrons in small balls of radius 8', R z s/4 as 
close as possible to each nucleus. We put the electron on top of the nucleus if the 
nucleus is not too close to the boundary of D \ D' and we place it closeby otherwise, 
thanks to the cone property. When the ball sits on top of the nucleus we call this 
a "perfectly screened nucleus" whereas we call the other ones "dipoles" . For later 
purposes we have to make sure that only the nuclei which are very close to the 
boundary are not completely screened. So we choose 

S' R z = min (<5 fljZ ,£/20). 

The factor 1/20 has no real significance but it is here to ensure that in any cone of 
size e — e' which is completely enclosed in D\ D', there is always a ball of radius 
e/5 in which there cannot be any dipole. 

Lastly, we have to cope with the fact that the charges in D' close to the boundary 
of D' do not have their optimal charge z opt , but rather the normal charge z = 
■2 op t + Sz. This additional positive charge 5z might create important electrostatic 
errors and we also have to screen it by adding electrons outside of D' . In spirit we 
follow the technique of 17 . To any unit cell Wj which is at a distance ^ e to the 
boundary of D' , we associate a little cone of size e in D \ D' , at a distance ^ e 
to Wj. This cone only depends on D and D', it does not depend on the random 
variable u. In this cone we know that there is a ball of radius e/5 at a distance 
> e/5 to the boundary of the cone, hence also at a distance ^ e/5 to the boundary 
of D 1 . We put the screening electrons in a small ball Bj of fixed radius ~ e in the 
cone, at a distance ^ e/5 to the boundary of the cone. Their total charge must be 
equal to 

5Zj := 5z - 

{R,z)eKC\W 3 

Note that each cone can intersect a finite (bounded) number of the other cones. 
This is because the cells Wj are at distance < 2e to their corresponding cone. So in 
the cone we can always reduce the size of the balls in which we put the electrons, 
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to make them all fit without any overlap. On the contrary to [17] where there was 
a smallest distance between all the nuclei, in our situation the additional electrons 
cannot always be chosen at a finite distance from all the other charges. However our 
construction guarantees that they can only overlap with perfectly screened nuclei, 
never with dipoles. 

If we summarize the situation, we have in our system (see Figure [3]) 

• electrons in D' , chosen to minimize the energy T(lo,D'), with the optimal 
charges close to the boundary; 

• nuclei in D' . They have a charge which might be larger than the optimal 
one when they are close to the boundary of D'; 

• classical dipoles outside of D' , at a finite distance ^ e/10 to the boundary 
of D\D'; 

• electrons in balls of a fixed radius, with a charge SZj used to compensate 
the charges of some of the nuclei in Wj C D' . They are at a distance ^ 2e 
but e/5 to the boundary of D' . They can never overlap with the dipoles; 

• perfectly screened nuclei living in D \ D', at a distance at least e/10 to 
the boundaries of D and D 1 . They do not interact with anybody, except 
possibly with the radial electron which we might have added in order to 
compensate some charges in D' . 





Figure 3. The trial state used in the proof of Lemma IT6l The 
picture is here for a cubic lattice with always exactly one nucleus 
per unit cell. 

(ii) Estimates on the energy of the trial state. We call J the set of all the indices 
such that Wj C\ (D \ D') ^ 0, dj the set of indices such that Wj intersects the 
boundary of D \ D' , and d Jo the set of indices such that Wj contains a nucleus 
whose charge has been optimized. For any j £ OJq, there is a ball Bj outside of 
D', containing the additional electrons with total charge SZj. 

Now we estimate all the terms. First, we have to pay for the kinetic energy to 
put the additional electrons in D\D' . Recall that we have two kinds of electrons, 
those which are squeezed in small balls of radii S' R and those which are used to 
compensate the charges of the nuclei close to the boundary of D 1 . The latter live 



THERMODYNAMIC LIMIT FOR DISORDERED QUANTUM COULOMB SYSTEMS 23 



in a ball of a fixed radius. The total kinetic energy of all these electrons, which we 
denote by K^ D ,, is bounded from above by 

jej (R,z)eicnWj v R > z> jej 

j&J j&Jo 

This inequality has been obtained by using the first eigenstates of — A on each ball, 
computing the associated Hartree-Fock state, and using it to define the electronic 
state in each ball. Note the power 5/3 which accounts for the fermionic nature of 
the electrons. Taking the average and using the regularity properties of D and D', 
we obtain the bound 

e ( k d\d>) < C\D\D'\ EI^ + C\D'\ 2 ' 3 El 5/3 < C(\D\D'\ + \D\ 2 ^ 3 ). 

Except for the kinetic energy, the nuclei which are completely screened do not 
participate much in our system, by Newton's theorem. There is only the possibility 
that they overlap with some electron in a ball Bj with j G J7b- Let us denote 
by pj the corresponding electronic density in the ball Bj and by Vr the Coulomb 
potential induced by the nucleus (R,z) together with its electron of size 5', Rz J^. 
We estimate the interaction between them using that each Vr vanishes outside of 
the ball of radius 5' R z /4: 



(R,z)eKnWk i 3 UR,z)eK.nB j nw k 




z 5 ' 2 

5/2 



X 



^CX^lo) 5 ' 3 + C Yl j 

(R,z)eicnB j nw k b ( s 'rJ/ a 

Ix ( Tj uj) 5 / 3 + X o(r k uj) 
\ w k nB^ 

Here we have used the Lieb-Thirring inequality to control J B p^ 3 by the kinetic 

energy, which in turn is bounded by Xo(rjUj) 5 ^ 3 . We get a similar term for every 
j € Jo- Summing over such j's, we deduce that the average of this error term is 
bounded above by a constant times \D\ 2 / 3 (E,X^ 3 +Elo). 

Our conclusion is that, up to an error of the form \D \ D'\ + |£>| 2 / 3 , we get the 
energy of a system in which we only have the quantum electrons in D' , interacting 
with classical particles. These are the nuclei in D' (with charges which might be 
higher than the optimal ones), as well as classical charges outside of D' . The latter 
are the dipoles of charges z and — z at a distance < e/10 to the boundaries of D and 
D', plus the additional electrons of charges SZj, used to compensate some charges 
in D'. We can write 

£(■,£>) < g+Z(-,D') + C\D\D'\ + C\D\ 2 / 3 , 
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with 

(47) Q=Z{;D')- f pr(W in + W out )+ z(W in (R) + W out (R))+IcUss, 

(B„z)eKnD' 

where W ou t is the Coulomb potential induced by all the classical particles sitting 
outside of D', and 

jeJo (R,z)eicnw 1 1 1 
is the potential corresponding to the excess charges of the nuclei in D' . Finally, 
/class is the Coulomb interaction between all these classical charges. 
In order to estimate Q in (j47|) . we use the method of pT]. We write 

) 



£(;D')-~ [ Pr (W in + W out 
V J 

- V z ( Wia (R) + W out (R))-^f 
1 (R,z)eKnD' 1 

VZ(;D')+ [1 + - ] /class- 



7 h 

The interpretation of the term in square bracket is that we have multiplied the 
charge of all the particles outside of D' by a factor —1/rj. Similarly we have changed 
the charges of the particles in D' to z opt — Sz/rj, instead of z = z opt + Sz. Now we 
claim that there is a stability of matter estimate in the form 



(48) E 



E(;D') + - [ pr(W in + W out ) 
V J 

- z(W in (R) + W out (R)) + - 



n £ — ' n 2 

(R,z)£KnD' ' 



\D\V S 



We first explain how to use this estimate, before turning to its proof in Step (in). 
Inserting this and using that £(•, D') ^ C\D\ by Lemma [15j we get an estimate on 
the average of Q: 

E(0) < E(£(-, D')) + C V \D\ + C 1 -^- + C E ( ^ 



rf/ 2 r] 

Using that the dipole-dipole interaction decays like R~ 3 at infinity, the classical 
interaction term can be estimated following the proof of Lemma IT51 (see Eq. (|43|0 : 

E(W) <C|£>| 2 / 3 log(|D|). 

The final estimate on (l47l) is 



m < Ha.*)) + am + clsgL + c w"^\ 

Now if we optimize in r/, finding that r\ oc |_D| -2 / 15 , and put back the other error 
terms, we arrive at our final estimate 

(49) E(Z(-,D)) < E(F(;D')) + C\D \ D'\ + C\D\ 13 / 15 . 

(Hi) Proof of the stability of matter estimate (j48|) . The estimate (|48|) is the equiva- 
lent of [17l Lemma 10] but, unfortunately, it does not follow from this result directly, 
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because there it was again assumed that there is a smaller distance between all the 
nuclei. To cope with this issue we slightly change the argument of [17j . 

The first step is the same as in [T7] and it consists in replacing the Coulomb 
potential l/|x| between all the particles by the Yukawa potential e - l x l/|x|. First 
we know that the Fourier transform of the difference is positive: 

1 1 



Ifcl 2 l + |fc| 2 



> 



and second we know that |a;| 1 — \x\ 1 e ' x ' — > 1 when \x\ — > 0. All this implies 
that for any yi £ M 3 and any charges ijjGl. 



E 



I 1 e-lw-f*, v „ 

v \Vi-Vj\ m-yj\ 1 / — 



2 



Thus, when we replace the Coulomb potential by the Yukawa potential, we get (in 
average) an error term of the form 

C[\D\- U 



2 

The first term is an estimate on the average number of electrons as well as the 
average number of nuclei in D' . The second term is an estimate on the average 
number of classical particles which are close to the boundaries of D and D', and 
whose charge has been multiplied by —1/?/. 

Now that we have replaced the Coulomb interaction by Yukawa, the second step 
of the proof consists in dropping all the classical negative charges. In a lower bound 
we only pay for the interaction with the positive charges. Consider for instance the 
interaction between the negative classical particles 1 jr\ inside or outside D' and the 
nuclei in D' which have the normal charge z or the optimized charge z opt . Because 
the Yukawa potential decays very fast, this interaction is easily controlled. It can 
be estimated similarly as in the proof of Lemma [15] by 

^ E E J- + r l~ 1 E J2 X ^)Xo(r k u)e-^- k \ 
j&(dj)(R,z)eWj R - z je(dj)kej 

almost surely. The first term accounts for nuclei which are in neighboring cells, 
whereas the second one deals with nuclei which are in non-adjacent cells. Hence the 
average of this term is bounded above by r?" 1 (E(Xi) +E(xf))#dJ < Cr]' 1 ^^ 3 . 
The argument is the same for the interaction between positive and negative charges 
r\~ x except that a crude bound gives Cr/- 2 \D\ 2 ^. All in all, we see that when we 
throw away the negative charges, we make an error which is bounded from below 
by -C|D| 2 /V 2 . 

Now we have reduced ourselves to a system of electrons interacting with nuclei 
through the Yukawa potential, up to a total error of the form -C(\D\ + \D\ 2 l 3 r 1 ~ 2 ). 
The nuclei in D' have a normal charge but the ones outside of D 1 have the charge 
z/rj. At this step we use the stability of matter with Yukawa, as was proved by 
Conlon, Lieb and Yau in 7 through the Thomas-Fermi Yukawa energy. First we use 
the Lieb-Thirring inequality and the Lieb-Oxford-type bound Eq. (A. 17)], and 
estimate the quantum energy from below by Thomas-Fermi theory. Then we use 
the lower bound on the form —CN — CJ2i z i^ 2 which is proved in [7, Eq. (A. 15)]. 
Hence, in average we get a lower bound of the form — C\D\ — C\D\ 2 / 3 r]~ 5 / 2 . Recall 
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that the other error terms are not worse than \D\ 2 / 3 rj 2 , which is itself smaller than 
1^12/3^-5/2 for v<1 

Now that we have proved (|48)) , this concludes the proof of Lemma [16j □ 

Remark 17. Using a Yukawa potential with mass fi and optimizing with respect 
to this mass in the end, it is possible to improve the error term |D| 13 / 15 . 

Step 3. Lower bounds. We already have one important lower bound on J-t,h 
and F_rp the one (|17|) corresponding to the stability of matter. This lower bound 
is true independently of to and of the shape of the domain D, which need not be 
regular. The existence of the thermodynamic limit (for simplices at least) follows 
from a much more precise lower bound which is stated in the following lemma. 

Lemma 18 (Graf-Schenker type inequality). Let A C R 3 be a fixed simplex (a 
tetrahedron). Then we have the following lower bound 

c\ f z t ^(u,d n g £A) 



(50) ^°)>^-j)] G -^ m ' *> 

c 



j(#{(R,z) eK(u)nD} + \D\) 



for every domain D, every t ^ 1, and with a universal constant C which only 
depends on the chosen simplex A. 

We recall that G — R 3 xi 5(9(3) is the group of translations and rotations acting 
on R 3 , endowed with its Haar measure. Since we have by convention T_t ^(ui, 0) = 0, 
the integral over G in the right side of (IBTfl) can be restricted to a compact set and 
it is therefore convergent. The proof of (1501) is based on an important inequality of 
Graf and Schenker 15; dealing with the Coulomb interaction of classical charges. 
This inequality was itself inspired of previous work by Conlon, Lieb and Yau [7J [H] 
and it is recalled in [17j Sec. 1.1.2]. 

The proof of Lemma ITSl is almost identical to that of [15] and [UJ p. 505-507]. 
We will not detail it again. The term z) g /C(w)nZ)} comes from the control 

of a localization error term in the Graf-Schenker inequality [15] . This estimate is 
conveniently done by using a version of stability of matter with Yukawa potentials 
which was derived in [Jj Eq. (A5)], based on a Thomas- Fermi-type theory. In [17] 
this error term was estimated by Z 5 / 2 times the number of nuclei, the latter being 
itself bounded by C\D\ (when D is regular). Here we do not have a smallest 
distance between the nuclei and we just keep the total number of nuclei in D. In 
this proof localizing the system to the small simplices gt A induces a small change 
of the nuclear charges close to the boundary of these simplices. This is why it is 
convenient to use the modified ener gy T_t . 

Remark 19. The error term #{(_R, z) £ K,{ui) PI D} is not necessarily bounded 
uniformly with respect to lo, but it has a limit almost- surely and in L P (R 3 ), when 
we divide it by the 'regularized' volume of D. Let us quickly explain this. For any 
a, e > 0, we introduce similarly to |16[ Eq. (9)] 



(51) 



|£>| a , e :=mf{|iy| : D' € U a , e , D'dD}. 



THERMODYNAMIC LIMIT FOR DISORDERED QUANTUM COULOMB SYSTEMS 27 



One can verify that this volume is of the same order as the volume of the union of 
the sets W + j with j € _£? intersecting D. Hence, for any fixed a, e > 0, we have 

#{(ii,z)e/CMn£>} < £ ]T 1= ]T Xo(r^) 

ie^f : (R,*)eK:(w)nW+7 ; 

(recall (|42p ). We /iawe A"o € L p (f2), according to assumption (|20|) and Lemma [7^] 
By i/ie ergodic theorem fTheorem \22\ below), this term behaves like \D\ a<s almost- 
surely and in L p (fl). In average we have the exact inequality 

(52) E ^ ^o(rjw) CE(Xq) \D\ a , e . 

(w+j)nD^% 

The inequality (j50|) (with a different error term) was denoted by (A5) in [TH] 
117] . There is a more precise but much more complicated inequality (A6) which is 
explained at length in these works. This inequality is also true in our case, provided 
we take the expectation value. 

Lemma 20 (A more precise lower bound for simplices) . Let D — gLA be a dilated, 
rotated and translated simplex. Then the property (A6) of [161 117] is valid for the 
averaged free energy E(.F T „(•, D)) . 

The proof is again exactly the same as in [17] . Instead of using the inequality [17] 
Eq. (89)] on the total number of electrons and on the kinetic energy in D, one 
uses 



From all the previous results we can now deduce that the thermodynamic limit 
exists for the deterministic function E(J[ T ) , by simply applying the abstract main 
theorem of [16]. 

Corollary 21 (Thermodynamic limit in average). There exists a function f(T,/j,) 
such that, for any sequence (D n ) C TZ a ,e of regular domains with a, e > 0, \D n \ — > oo 
and diam(D„)|Z?„| -1 / 3 ^ C, we have 

Efj^T J- ,D n )) 
(53) lim 1 in = /(?»• 

n^oo \L) n \ 



Proof. The function D H> E(J 7 t. a ,( • , -D)) satisfies all the properties (A1)-(A6) 
of Q£]. For (A3) it is even periodic, E(J" T , M ( ■ ,D + k)) — E(j r TjAI ( • , £>)), Vfc e _Sf, 
by stationarity. We deduce the result from [Jj)J Thm. 2]. □ 

Our task in the next step will be to upgrade the convergence of the average to 
a convergence in L q (fl). Since we already know from ([571) that \D\~ l T_ T „(•, D) is 
bounded in D>/ 2 {£1), strong convergence in L l (£l) will imply the result by interpo- 
lation. 

Step 4. Thermodynamic limit in L 1 (il). This step is now more specific to the 
stochastic case. We will show that for a sequence (D n ) of regular domains like in 
the statement of Theorem^ \Dn\~ 1 F_ T ^ , D n ) converges strongly in L^fi) to the 
same limit as its average, namely f(T,/j,). 

We will make use of the following celebrated theorem. 
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Theorem 22 (Ergodic theorem [H]). Let X be a random variable in L p (£l) for 
some 1 ^ p < oo. Consider a sequence of sets D n C R 3 such that \D n \ — > oo and 
which is regular in the sense of Fisher [13] . that is 

(54) Vte[0,to], \{xeR 3 : d(x,dD n ) <t|£> n | 1/3 }| < \D n \n{t) 

for some to > and some function r\ with lim^o vW) — 0- Then we have 

1 



lim E 



= 0. 



fce.S?n.D„ 

If moreover D n C B c i£) 11/3 /or some c > and aZ/ n, then 

1 



^ X{ Tk u) — > E(X) 



almost- surely. 



Theorem [22] can for instance be found in a paper of Tempel'man [34]. The 
Fisher regularity assumption is stronger than Tempel'man's condition (E\). The 
convergence in L p (f2) is Theorem 6.4' of [34], whereas the almost-sure convergence 
is Theorem 6.1 in [34] . 

Our first useful result is the following 

Lemma 23 (Limit of negative part). Let (D n ) C 7Z a ,e be a sequence of regular 
domains like in the statement of Theorem^ Then we have 



(55) 



lim E 



\D r , 



-f(T,n) 



= 0. 



Proof. The proof mainly follows from the Graf-Schenker inequality (|50p which is 
our only precise bound valid almost-surely. We know that for all i 1 



T T Juj,D n ) > l-C/£ f Z T>IM (uj,D n n(g£A)) 



\Dn\ 



IAJ 



\£A\ 
C 



■dg 



D n 



■#{(R,z)€lC(u)nD n } 



Note that we can write, as mentioned in Remark [19] and using the regularity of D ni 

1 #{(fi,2)e/CHno4^^ioM 



\D 



where J n is the set of all points j of the lattice ££ such that (W+j) n D n ^ 0. By 
the ergodic theorem (Theorem |2"2"]) and Lemma IT41 we know that 



(56) 



lim E 

n— >oc 



J- X o(Tju) E(X ) 

jEJn 



0. 



The limit also holds almost-surely provided that D n C B c \ D n/3. 
particular that 

£r, M (w,D„) > 1 - C/£ f ZtA u > d » n ^)) ^ C 



This means in 



I A, 



|*A| 



4ff-y(l + y»), 
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where Y n — > strongly in L 1 (SI). Now we argue similarly as in p~6j Proof of Thm. 1]. 
We write the integral over translations and rotations as 



(57) 



1 r Z. T ^(u,,D n n(g£A)) 



\D 



dg 



n I JG 



\£A 



D n 



dR 



W JSO(3) 



Z T>ll (uj,D n n(MA + T + j)) 
\£A\ 



We denote by Q n the set of all j G if such that R£A +t + j C D n for all t e W 
and all R € SO (3). Similarly we denote by dQ n the set of all the indices such that 
(R£A + r) + j n dD n ^ for at least one i? e 50(3) or one r e W. Using the 
regularity of D n we have #dQ n ^ C|Z3 n | 2 / 3 and we obtain that 



\D r . 



\D n ng£A\ 
\£A\ 



dg 



= \W\(#Gn)+ V / dr f dR 

:^a,n JW JSO(3) 



jedg 



D n n(R£A + T + j)\ 
\£A\ 



^ \w\(#g n ) + \w\(#dg n ) ^ \w\(#g n ) + c\D n \ 2 '\ 

see [H Eq. (17)]. Hence #£„ = (W^l-Onl + o(\D n \). Coming back to ([57]) and 
using the stability of matter 

Z T>ll (u,D n n {R£A + r + j)) > -C\£A\ 

for the indices j £ dG n , we get 



(58) 



1 r Z T Joj,D n n(g£A)) 



D 



71 1 JG 



\w\ 



dr 

W JSO(3) 



\£A\ 
dR 



dg 



l + o(l) y, Z^fow.itfA + T)) 



*Gn 



jeQ„ 



\£A\ 



Therefore we have proved that 
Z T Ju,D n ) > l-C/£ 



\D„ 



\W\ 



dr 



dR\ 



W JSO{3) 



l + o(l) y, T t Jt j u 1 R£A + t)) 



\£A\ 
C 



C\D n \- 1/3 -j(l + Y n ). 



We now subtract f(T, /i) on both sides, take the negative part and average over uj. 
We obtain 



E 



£-T,tS.'t Dn) 



C 



< C|D n |-V3 + i x +nYnl +/(T)/i) ) 



1 



-(l-C/£)— - / dr 



dRE 



W JSO(3) 



i + o(i) ^2 £ T A T J-> mA + T )) _ /(T M) 



je6?„ 



|*A| 



30 



X. BLANC AND M. LEWIN 



By the ergodic theorem (Theorem [22]) , we have 



lim 



l + o(l) £r,„fc-.-RlA + T)) _ E J r J-,itfA + r)) 



jeg„ 



\£A\ 



\£A\ 



By the dominated convergence theorem, we deduce that 



lim -— / dr / dRE 

n^oo \ W\ J w J SO (3) 



1 

W\ 



1 + (1) ^ Zt, m (^-.^A + t)) 



IV J SO{Z) 



|£A| 

7" Tu (-,i?IA + T)) 

e= |7a1 ~-^ T ^) 



Thus, 

(59) lim sup E 

n— >oo 

l 



I A, I 



€ -. — r I dr dR 

\W\ J w 7 SO (3) 



E 



£r, ,(^A+r)) 
|£A| 



+ ^(l + /(T,/x)). 



We have proved the existence of the thermodynamic limit for EF_ T „ in Coroilarv[2"Tl 
On the other hand we have a uniform bound 



E 



by Lemma [T51 By the dominated convergence theorem, we conclude that the right 
side of (1551) tends to when I — > oo , hence that 



lim E 

n—>oo 



l-DJ 



This concludes the proof of Lemma [ 



□ 



Remark 24. If we assume that D n C B c \ Dn y./3, then we can actually show that 

J-rp D n ) 

(60) lim inf ' ' >f(T,fx) 



Dr. 



almost surely. Indeed, we have 



jeg„ 



\£A\ 
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by the stability of matter (Theorem Coming back to ()58[) . we can now use 
Fatou's Lemma and the almost-sure ergodic theorem (Theorem \22\) to infer 



lim inf • 



1 f Z Tjfl (u>,D n n(g£A)) 



D 



n\ Jg 
1 



\w\ 
1 

W\ 



dr 



w 



SO(3) 



dg 



dR lim inf 

n— >co 



1 + 0(1) ^2 ^( T ^: ReA + T )) 



dr 



w 



dR E 



SO(3) 



\£A\ 



\£A\ 



The proof of (|60l) is then the same as before. 

Now we can deduce our desired result. 

Corollary 25. Let (D n ) C TZ a .e be a sequence of regular domains like in the 
statement of Theorem^ Then we have 



(61) 



lim E 

n— >oo 



\D„ 



= 



for q — 1 if p — 2 and all 1 ^ q < p/2 if p > 2. 



Proof. We write 

£t,/j('i^i) 



E 



Ad 



sEt,j*(*> Ai) 



= E: 



I Ail 



/(2»+2I 



Zt,u(-,A) 



The term on the right converges to by Corollary [21] and Lemma [23l Thus we get 
convergence in L l (Vl). In addition, we know that Z-Tai'i D n )\D n \~ l is bounded 
in L p / 2 {Q) by Lemma [T5l Hence, by interpolation, we have convergence in L q (fl) 
with q = 1 if p = 2 and Kg< p/2 if p > 2. □ 



Remark 26. If D n C Z? c i£) 11/3, i/ien we actually deduce that 
(62) lim inf = /(T, /x) 



n->oo |D„| 

Indeed, we have by Fatou's Lemma (using that there is a uniform lower bound by 
the stability of matter) 



.Zt, m ( w > A) 



A 



sC E I lim inf ^^l', - - f(T,fi) ) < liminf ( E ^'*""' "" - f(T,n) ) = 



,Zt,^( w > A) 



An 



and (|62p follows. 
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Step 5. Back to the original free energy. Our last step is to come back to 
the original energy Tt,^{<^->, D) for which we do not optimize over the charges of the 
nuclei which are close to the boundary. First we have J-y !jU (o;, D n ) ^ J T(1 (u,I) n ) 
by definition, hence obviously 



^T,fi{-,D n ) 



I 



f{T,n) 



^ E 



A 



0. 



by Lemma [221 Consider now a sequence (D' n ) C TZ a ', E ' such that D' n C D n , 
\ D n\ = IA| + o(\D n \) and d{dD' nl dD n ) ^ 4e. The sets D' n can for instance be 
constructed by considering a tiling of K 3 and taking the union of all the cells in 
D n which are at the appropriate distance to the boundary of D n . The regularity 
of such sets is proved in [16| Prop. 2]. By Lemma [TBI we have 



E 



\D rj 



< E 



which proves that 



lim E 

n— >oo 



\ |A| 
\ |A| 



s=E 



V \D' r , 



o(l) 



The convergence in L q (O) follows from the proof of Corollary [ 
the proof of Theorem 0J 

5. Proof of Corollary [5j 



This concludes 
□ 



Let us denote by T n (w) the optimal electronic state in Fock space for Tt^ (w, D n ) , 
and by I n (u) the total Coulomb interaction energy between the nuclei and the elec- 
trons in D n . These quantities depend on T and p but we do not emphasize this 
in our notation. We have by the variational characterization (|16p of J~t,h and by 
Lemma [15] 



^J„ + Tr ( ^(-A)j - pM ) r„ + T Tr^ T n log r„^J = E F T . t 



{;D n )^C\D n \. 



We know from (XTl Lemma 9] that 



Tr r„ +TTr^r n lo g r„ > -Tlog (W e"( E,(-A),-2/^)/ 

> -C\D n \. 

By the Lieb-Oxford inequality [351 OSO and the estimate on E J p^/ 3 + pr n , we 
have 

•: Tr ( Y - 



r„ > -ED(p r J 



p^ 3 > -ED(prJ-C|D n | 



with 



D(pr 



Pr„(x)pr„(?/) 



'D„jD n \x-y\ 

denoting the classical Coulomb energy. We deduce that 
'1 



dx dy 



E 



V nP T n +U n ) ^C\D n \ 



D„ 
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where V n {uj, x) is the electrostatic potential induced by the nuclei in D n and U n {uj) 
is their Coulomb repulsion. At this step we replace each nucleus (R, z) in D n nlC(uj) 
by a smooth, spherically symmetric, distribution vr leaving in a ball of radius 
min(l, 5r iZ /3), the distance to the closest other nucleus. By Newton's theorem, U n 
does not change and the electrostatic potential V n (uj,x) is modified only when x is 
in these balls. The new potential V n (u>,x) satisfies 



V n (u,x) - V n {uj,x) 



<* E 

(R,z)£D n nK(u>) 



Using that 



PvAVn-Vn 







5/2 


3 r 




V n - V„ 


H 




5 JD n 




5Jd 



5/3 



and the estimate 

V n ~ Vn 



5/2 <z 5/2 



E 



(i?, Z )GD„n/C(a;) JS ( i? ' min ( 1 ' ,5 «-/ 3 )) 

<s cz 5/2 E L 

(R,z)eD n nK.(uj) 



dx 
|x| 5 / 2 



we deduce from the Lieb-Thirring bound ([38)) and the fact that EJq < oo, that 



2 E D I Pr„ 



]T zv R \ ^c\d„\ + ~eI 

(R,z)£D n nK(-) J \ 



E ^°(^) 

(i?,z)en„n/c 



Since -D (vr) — C/Sr <z and E(Xi) < oo by assumption, the last term is bounded 
by a constant times \D n \ and we end up with the bound 



\ {R,z)£D n nlC(-) 

At this step we use a capacity estimate, namely 



D(/) > ( [ A inf D(g) >( [ f) inf D 



(9) = CT: 



diam(Z)„) 



fa=i 



where B n is the smallest ball containing D„ + B\. Using the regularity of D n as 
well as the fact that diam(D„) ^ C| Z?^ | 17,3 , we get our final estimate 



E 



\D 



n\ JD n 



E 



(R,z)£D n nK(-) 



C 



|A>| 2/3 ' 



Since Xq £ L 2 (f2) and z £ L°°(f2) by assumption, we have by the ergodic theorem 
(Theorem W$ 

2 

1 



lim E 

n— >oo 



E 



z - Z» 



(R,z)GD n nK(-) 



= 
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and it follows that 



lim E 



1 



= 0. 



I Ail JD n 

This is valid for any chosen temperature T ^ and chemical potential /i. 
To conclude the proof, we use the variational principle (fT())) and get 

Jr,o(w,Ai) - ^N Tttl (u>, D n ) < TT,n{u,D n ) ^ T T ,o(uj,D n ) - p N Tj0 (u), D n ), 

where we recall that Nt^{u>, Ai) — / Pr„ is the number of electrons at temperature 
T and chemical potential /i. From the limit of Nt,u(', > AO m ^ 2 (^)j hence in 
i 1 (r2), we finally conclude that 

,-?T,/i(-, Ai) 



lim E- 



I A 



,. m -^T,o(-! Ai) „ 

hm E — — /iZ a 

n-s-oo L)„ 



f(T,0)-fiZ a 



It is clear from the variational principle that T t— > E 7r,o('i-0n) is concave, hence 
the limit /(T, 0) must also be concave. This concludes the proof of Corollary[5j □ 

Let us end this paper by two final remarks. 

Remark 27 (A better convergence for simplices). Like in [T6l |T7] 7 it is possible 
to strengthen the result when D n is a sequence made of dilated simplices which are 
possibly rotated and translated. More precisely, we have 



lim 



\£A\ 



f{T,ii) 







L*(n,L°°(G)) 



for all 1 ^ q < p/2 if p > 2 and q = 1 if p = 2. 



Remark 28 (An abstract setting). Our approach in this paper is general and it can 
be stated in an abstract fashion, like in |16j . Consider a random variable J-(u),D) 
defined on all bounded open subsets D ofW 3 . Assume that J- satisfies the following 
assumptions: 

(AI) F(u), 0) = G a.s.; 

(A2) J~(uj, D) — C for all D and almost all to; 

(A3) J~(w, D + k) = J-(TkU>, D) for all D, all k € ££ and almost all lu; 
1 - a(l) 



(A5) F(u, D) > - ~ y ~' I D n glA) dg - (\D\ a , e + R D {ia))a{l) for all 

D and almost all lu, where Rd is such that E|_Rd| ^ C|Z?| Q . £ for all D; 

(A4-6) D i— > EJ r (-,D) satisfies the assumptions (A4) and (A6) of [16] . 

Under these conditions, the results of j!6j as well as the method of the present 
article show that there exists a constant f such that 

J"(-,A 



lim 



I A I 



-/ 



= 



for any growing regular sequence (D n ) like in Theorem^ Also 
when D„ C B c i D n/ 3 . 



. -F(avA) 
hmml — — — | — = f a.s. 

n-¥oo \L> n \ 



THERMODYNAMIC LIMIT FOR DISORDERED QUANTUM COULOMB SYSTEMS 35 



References 

[1] M. AlZENMAN AND S. Warzel, Localization bounds for multiparticle systems, Comm. Math. 

Phys., 290 (2009), pp. 903-934. 
[2] P. W. Anderson, Absence of diffusion in certain random lattices, Phys. Rev., 109 (1958), 

pp. 1492-1505. 

[3] D. Basko, I. Aleiner, and B. Altshuler, Metal-insulator transition in a weakly interacting 
many-electron system with localized single- particle states, Annals of Physics, 321 (2006), 
pp. 1126 - 1205. 

[4] J. R. Baxter, Inequalities for potentials of particle systems, Illinois J. Math., 24 (1980), 
pp. 645-652. 

[5] X. Blanc, C. L. Bris, and P.-L. Lions, The energy of some microscopic stochastic lattices, 
Arch. Ration. Mcch. Anal., 184 (2007), pp. 303-339. 

[6] V. Chulaevsky and Y. Suhov, Multi-particle Anderson localisation: induction on the num- 
ber of particles, Math. Phys. Anal. Gcom., 12 (2009), pp. 117-139. 

[7] J. G. CONLON, E. H. Lieb, and H.-T. Yau, The N 7 / 5 law for charged bosons, Commun. 
Math. Phys., 116 (1988), pp. 417-448. 

[8] , The Coulomb gas at low temperature and low density, Commun. Math. Phys., 125 

(1989), pp. 153-180. 

[9] B. Deissler, M. Zaccanti, G. Roati, C. D'Errico, M. Fattori, M. Modugno, G. Mod- 
UGNO, AND M. INGUSCIO, Derealization of a disordered bosonic system by repulsive interac- 
tions, Nature Physics, 6 (2010), pp. 354-358. 
[10] F. J. Dyson, Ground-state energy of a finite system of charged particles, J. Math. Phys., 8 
(1967), pp. 1538-1545. 

[11] F. J. Dyson and A. Lenard, Stability of matter. I, J. Math. Phys., 8 (1967), pp. 423-434. 
[12] C. Fefferman, The thermodynamic limit for a crystal, Commun. Math. Phys., 98 (1985), 
pp. 289-311. 

[13] M. E. Fisher, The free energy of a macroscopic system, Arch. Ration. Mech. Anal., 17 
(1964), pp. 377-410. 

[14] M. E. Fisher, and D. Ruelle, The stability of many-particle systems, J. Math. Phys., 7 
(1966), pp. 260-270. 

[15] G. M. Graf and D. Schenker, On the molecular limit of Coulomb gases, Commun. Math. 
Phys., 174 (1995), pp. 215-227. 

[16] C. Hainzl, M. Lewin, and J. P. Solovej, The thermodynamic limit of quantum Coulomb 
systems. Part I. General theory, Advances in Math., 221 (2009), pp. 454-487. 

[17] , The thermodynamic limit of quantum Coulomb systems. Part II. Applications, Ad- 
vances in Math., 221 (2009), pp. 488-546. 

[18] P. D. HlSLOP, Lectures on random Schrodinger operators, in Fourth Summer School in Anal- 
ysis and Mathematical Physics, vol. 476 of Contemp. Math., Amer. Math. Soc, Providence, 
RI, 2008, pp. 41-131. 

[19] W. KlRSCH, An invitation to random Schrodinger operators, in Random Schrodinger op- 
erators, vol. 25 of Panor. Syntheses, Soc. Math. France, Paris, 2008, pp. 1-119. With an 
appendix by Frederic Klopp. 

[20] F. Klopp and B. Metzger, The Gross- Pitaevskii Functional with a Random Background 
Potential and Condensation in the Single Particle Ground State, Math. Phys. Anal. Geom., 
pp. 1-21. 10.1007/sll040-011-9099-x. 

[21] A. Lenard and F. J. Dyson, Stability of matter. II, J. Math. Phys., 9 (1968), pp. 698-711. 

[22] E. H. Lieb, The stability of matter, Rev. Mod. Phys., 48 (1976), pp. 553-569. 

[23] , The stability of matter: from atoms to stars, Bull. Amer. Math. Soc. (N.S.), 22 (1990), 

pp. 1-49. 

[24] E. H. Lieb and J. L. Lebowitz, The constitution of matter: Existence of thermodynamics 
for systems composed of electrons and nuclei, Advances in Math., 9 (1972), pp. 316-398. 

[25] E. H. Lieb and S. Oxford, Improved lower bound on the indirect Coulomb energy, Int. J. 
Quantum Chem., 19 (1980), pp. 427-439. 

[26] E. H. Lieb and R. Seiringer, The Stability of Matter in Quantum Mechanics, Cambridge 
Univ. Press, 2010. 

[27] E. H. Lieb and W. E. Thirring, Bound on kinetic energy of fermions which proves stability 
of matter, Phys. Rev. Lett., 35 (1975), pp. 687-689. 



3(5 



X. BLANC AND M. LEWIN 



[28] E. H. Lieb AND H.-T. Yau, The stability and instability of relativistic matter, Commun. 

Math. Phys., 118 (1988), pp. 177-213. 
[29] L. Onsager, Electrostatic interaction of molecules, J. Phys. Chem., 43 (1939), pp. 189-196. 
[30] D. Ruelle, Classical statistical mechanics of a system of particles, Helv. Phys. Acta, 36 

(1963), pp. 183-197. 

[31] D. Ruelle, Statistical mechanics of quantum systems of particles, Helv. Phys. Acta, 36 
(1963), pp. 789-799. 

[32] D. Ruelle, Statistical mechanics. Rigorous results, Singapore: World Scientific. London: 

Imperial College Press , 1999. 
[33] L. Sanchez-Palencia, Quantum gases: Joint forces against disorder, Nature Physics, 6 

(2010), pp. 328-329. 

[34] A. A. Tempel'man, Ergodic theorems for general dynamical systems, Trudy Moskov. Mat. 

Obsc., 26 (1972), pp. 95-132. 
[35] N. A. Veniaminov, The Existence of the Thermodynamic Limit for the System of Interacting 

Quantum Particles in Random Media, arXiv:1112.2575 (2011). 

CEA, DAM, DIF, F-91297 Arpajon, France. 

E-mail address: blancOann.jussieu.fr, xavier.blancOcea.fr 



CNRS & Laboratoire de Mathematiques (UMR 8088), Universite de Cergy-Pontoise, 
F-95000 Cergy-Pontoise, France. 

E-mail address: mathieu.lewinOmath.cnrs.fr 



